Electron transport through an interacting region: The case of a nonorthogonal 

basis set 
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The formula derived by Meir and Wingreen [Phys. Rev. Lett. 68, 2512 (1992)] for the electron 
current through a confined, central region containing interactions is generalized to the case of a 
nonorthogonal basis set. As the original work, the present derivation is based on the nonequilibrium 
Keldysh formalism. By replacing the basis functions of the central region by the corresponding 
elements of the dual basis the lead- and central region-subspaces become mutually orthogonal. The 
current formula is then derived in this new basis, using a generalized version of second quantization 
and Green's function theory to handle the nonorthogonality within each of the regions. Finally, the 
appropriate nonorthogonal form of the perturbation series for the Green's function is established 
for the case of electron-electron and electron-phonon interactions in the central region. 

PACS numbers: 73.23.-b,73.63.-b,73.63.Nm 



I. INTRODUCTION 

Electron transport in nano-scale contacts is a highly 
active research area. During the last decade it 
has become possible to create two-terminal junctions 
where atomic-sized conductors are contacted by macro- 
scopic metal electrodes using scanning tunnelling mi- 
croscopes^, mechanically controlled break junctions^, 
or electromigration techniques 5 . In this way I-V char- 
acteristics have been obtained for a variety of differ- 
ent nano-contacts including carbon nanotubes 6 , metal- 
lic point contacts^, atomic wires^, as well as individual 
molecules ranging from large organic compounds^ down 
to a single hydrogen molecule^. 

The quantitative modelling of the electrical properties 
of a nano-scale contact represents a great theoretical chal- 
lenge involving a detailed description of both the atomic- 
arid electronic structure of a current-carrying system out 
of equilibrium. At present, the most popular approach 
to the problem combines a non-equilibrium Green's func- 
tion (NEGF) formalism with ab initio electronic struc- 
ture theory. A cornerstone of this approach is a formula 
giving the current through the system in terms of the 
Green's function of a spatial region containing the con- 
tact (the central region). When interactions are limited 
to the central region, the current formula is an exact 
result, however, in practice some approximation for the 
full interacting Green's function must be invoked. For 
example, in the commonly used NEGF-DFT approxima- 
tion the exact Green's function is replaced by the non- 
interacting Kohn-Sham Green's function defined within 
Density Functional Theo^ -A^U ,, 

Application of the NEGF theory to electronic trans- 
port in quantum wires was introduced as an alternative 
to the Landauer-Biittiker formalism to treat electronic 
interactions. In 1991 Hershfield and co- workers^ derived 
an expression for the current through a single interact- 
ing level (an Anderson impurity) , and the following year 
Meir and Wingreen^ generalized the current formula to 
the case of an arbitrary number of states in the central 



region [Eq. (6) of Ref . Il3| . In these studies the system, 
i.e. the conductor, was partitioned into three parts (the 
leads plus the central region) , and the basis of the single- 
particle Hilbcrt space was taken as an orthonormal set of 
functions each belonging to one of the three regions. In 
practical ab initio calculations, however, the requirement 
of localized basis functions, which is essential for the par- 
titioning, is difficult to combine with orthogonalitj*i2ii£. 
It is therefore of great practical importance to generalize 
the current formula to nonorthogonal basis sets. 

In this paper, a rigorous operational framework for ap- 
plying second quantization and Green's function methods 
in a nonorthogonal basis is presented and used derive a 
generalized current formula which is valid for nonorthog- 
onal basis sets. The main problem in deriving the cur- 
rent formula in the general case is the lack of orthogo- 
nality between basis functions in the central region and 
basis functions in the leads (due to the assumption of 
localized basis functions, we can take basis functions be- 
longing to different leads to be non-overlapping and thus 
orthogonal). We solve this problem by replacing the ba- 
sis functions of the central region by the so-called dual 
basis functions which, by construction, are orthogonal to 
the basis functions of the leads. The nonorthogonality 
of the basis within each of the three regions is handled 
by applying a generalized version of second quantization 
and Green's function theory. 

For non-interacting electrons, the NEGF current for- 
mula takes a particularly simple form [Eq. (7) of Ref. fl3| . 
This version of the formula, which is equivalent to the 
well-known Landauer-Biittiker formula, has been applied 
extensively over the last years and forms the basis of 
most numerical schemes addressing phase-coherent trans- 
port in nano-scale structures. In the special case of no 
interactions, Xue et al. have derived the form of the 
current formula in a nonorthogonal basis^lS Xue et al. 
takes a route which is less direct than the one by Meir 
and Wingreen as it involves transformations between the 
basis and real space representations. In contrast, the 
derivation given here follows the original work closely 
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and is formulated entirely within the basis set represen- 
tation. More importantly, the derivation of Xuc et at 
cannot be extended to the interacting case since it re- 
lies on a spectral representation of the Green's function 
in terms of single-particle eigenstates. Emberly and Kir- 
czenowi& have addressed the problem of nonorthogonal- 
ity within the Landauer-Buttiker formalism. Here the 
current due to non-interacting electrons is obtained di- 
rectly from the singe-particle scattering states which are 
evaluated in a Hilbert space with an energy-dependent 
inner product. Finally, Fransson et al. have studied 
the effects of nonorthogonality within the tunnelling for- 
malismiALiS Using the Kadanoff-Baym approach and the 
Hubbard operator technique, they derived an expression 
for the current through a single interacting level (Ander- 
son model) in the weak coupling limit, taking the finite 
overlap of the tunnelling wave functions into account ^ 
In contrast, the derivation presented here makes no as- 
sumptions about the type of interactions and is not lim- 
ited to the case of weak coupling between the leads and 
the central region. Moreover, the use of the dual basis in 
the central region simplifies the formalism significantly 
and avoids many technical problems otherwise arising 
from the nonorthogonality between the lead- and central 
region-states. 

The paper is organized as follows. In Scc.[H]thc general 
theory of second quantization and Green's functions in a 
nonorthogonal basis is reviewed, and the concept of the 
dual basis is introduced. In Sec. lIIII these results are used 
to derive the current formula in a nonorthogonal basis, 
both in the interacting and the non-interacting case. 



II. NONORTHOGONAL BASIS SETS 

In this section we generalize the second quantization 
formalism and elements of the Green's function theory 
to the case of nonorthogonal basis sets. To avoid mathe- 
matical problems with convergence, we assume that our 
single-particle Hilbert space, TL, is finite dimensional. 
While this assumption might seem unsatisfactory from 
a fundamental point of view, it will always be true in 
practical applications. 



A. Second quantization of one- and two-body 
operators 

Throughout this section let {&} denote a - not neces- 
sarily orthonormal - basis of the single-particle Hilbert 
space, TL. The corresponding creation and annihilation 
operators^, ej and Cj , acting on the fermionic Fock space 
fulfill the canonical anti-commutation relations^ 

[c\,c]} = (1) 
{c i)Ci } = (2) 
{ c ii c ]} — Sij, (3) 



where Sij = (<pi\<pj) is the overlap matrix. 

To any one-body operator, A^ x \ acting on TL, we asso- 
ciate the matrix Aij = (4>i\A^ \<f>j). We have the follow- 
ing representation of A^ in terms of the basis vectors 

{<M, 

Aw^sy&XM ( 4 ) 

ij 

where we have used Dirac's bra-ket notation and intro- 
duced the matrix 21 = S~ 1 AS~ 1 . This matrix trans- 
form will appear often throughout the text, and we shall 
reserve the Gothic font for matrices that results from 
such a transformation. The validity of the representa- 
tion is easily checked by evaluating the inner products 
{(f) n \A^ 1 '\(f) m ) on both sides of the equation. The second 
quantized form of A^\ which we denote by A, is 

y 

The easiest way to derive this expression is to start from 
the well known form of A in terms of some orthonormal 
basis, {ipn}-, and corresponding creation and annihilation 
operators d) n , d n and then expand these in terms of the 
original c|, Cj. The explicit form of this expansion reads 
V^n = J2ij{S~ 1 )ij{ ( l>j\' l Pn)\<f>i}, which follows from Eq. (gj) 
applied to the identity operator. 

We now turn to a general two-body operator, B^ 2 \ de- 
fined on the two-particle Hilbert space TL^ — TL ®TL. 
A basis for TL^ is provided by the tensor products 
{4>i ® <pj}. We define the corresponding overlap matrix 

Sij'kl = ® <t>j\4>k ® 4>i) = SikSji, as well as the ma- 
trix Bij t ki — {<t>i ® 4>j\B^ \4>k <S> 4>i)- As for the one-body 
operators we have the first quantized representation 

B m = By,/w|&®^)(fo<8>&|, (6) 

i,j,k,l 

where the matrix 23 is defined by S = 
(S^y 1 B(S ( ^)' 1 . For later use we note that 

(S (2) )£ M = (S- 1 ) ik (S- 1 ) jl , (7) 

which can be directly verified by multiplication with . 
Finally, the following expression for the second quantized 
version of can be derived using the same technique 
as for the one-body operator, 

B= %>i 3 .kic\c]cic k . (8) 



B. Single-particle Green's functions 

In order to fix notation we start with some well known 
definitions. Given two single-particle orbitals <pi and <pj 
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(not necessarily normalized or orthogonal) the retarded 
and advanced single-particle Green's functions (GFs) are 
defined by 

Gfc(t,i / ) = -i6(t-t')({ Ci (t),c](t')}) (9) 
G^(t,t') = l e(t'-t)({c l (t),c](t')}). (10) 

Here the brackets () denotes an expectation value with re- 
spect to the equilibrium state of the system. The greater 
and lesser GFs are defined by 

G<(M') = i(c){t') Ci {t)) (11) 

G>(t,t>) = -.^(iJcJCf)). (12) 

In fact all of these GFs can be derived from the contour- 
ordered GF which is defined by 

G ij (r,T') = -i{T c [c i (r)4(r')]). (13) 

Here t = (t, a) is a collection of the time variable t and a 
branch index, u, and Tq is the contour-ordering operator. 
Note that Cj(r) = Cj(i) and cj(r') = c](t'), while the 
branch indices merely serve to determine the ordering of 
the operators. For more comprehensive introductions to 
the general GF theory the reader is referred to Refs. EH 

We consider first the case where both the expectation 
value and the time-evolution of the creation and annihi- 
lation operators entering the GF is governed by a time- 
independent, quadratic Hamiltonian, h. Expressing h 
as in Eq. ©, using that dtCi(t) = i[h,Ci](t), and Fourier 
transforming with respect to the time difference t — t' , 
we obtain the equation of motion for the retarded GF 
matrix in the basis {4>i}'- 

(5" 1 o; + - S^hS-^Gr^) = I. (14) 

Here hij = (<j)i\h^ \4>j) where is the first-quantized 
version of h, and u) + = (u> + irj) with 77 a small positive 
number ensuring proper convergence of the Fourier inte- 
gral. The same equation holds for G a (u) when 77 — > —rj. 
It is useful to introduce another matrix quantity related 
to the GF and defined by 

<5 X =S- 1 G X S~ 1 . (15) 

As indicated by the superscript x the definition applies 
to any of the GFs introduced above. To have a name we 
shall refer to G5 as the overlap GF. Its retarded variant 
clearly fulfills the following matrix equation 

(Suj + -h)<8 r (u>) =1. (16) 

Next, we ask about the form of the perturbation series 
of the GF in a nonorthogonal basis. We thus consider a 
quantum system with a Hamiltonian H = h + V, where 
ho is a simple quadratic Hamiltonian, while V is a com- 
plicated one- or two-body perturbation. In the perturba- 
tion expansion of the contour-ordered Green's function, 



Gij(r, t'), we encounter the usual terms (generating dia- 
grams with two external vertices): 

J dri • ■■dr n (Tc[c i , ho (r)cl h y)V ho (r 1 ) ■ ■ ■ 14 (r n )]) . 

In this expression both the average and 
the time evolution is governed by h$, i.e. 
^o( T ) — enp(itho)X exp(— itho) and (^Qo = 
Tr[Xex P (-/?/\))]/Tr[exp(-^/in)], with = 1/kT. 
When V is represented in terms of a nonorthogonal 
basis as in Eqs. or (jHJ, the term (|17|) will generate 
unperturbed (n + l)-particle GFs (or (2n + l)-particle 
GFs if V is an interaction) involving creation and anni- 
hilation operators of the nonorthogonal orbitals <pi. As 
usual this (n + 1)- or [2n + l)-particle GF can be broken 
down into unperturbed single-particle GFs using Wick's 
theorem2ii22. Although Wick's theorem is normally 
proved for orthonormal states, its validity for nonorthog- 
onal states can be readily verified by expanding each 
creation/annihilation operator entering the (n + 1)- or 
(2n + l)-particle GF in terms of a fixed orthonormal 
basis. Wick's theorem can then be applied to each term 
in this expansion, which now involves only orthonormal 
states, and finally the original basis functions can be 
reintroduced. The perturbation series for the matrix 
Gij(r, t') in terms of the nonorthogonal basis, 
should therefore be constructed according the usual 
Feynman rules, using G° (t, t') as free propagator and 
?Gij (or 2Jy,fc() as the coupling strengths entering at 
the vertices. Equivalcntly, the perturbation series for 
the overlap GF, <Sij(r, t'), is obtained by evaluating 
each diagram using 6^(r, r') as free propagator and Vij 

(or Vij.ki) as coupling. In the case where V contains 
interactions this follows from the identity Eq. {7J . 

C. Dual basis 

Below we introduce the concept of a dual basis. The 
dual basis will be used in the next section for orthogo- 
nalizing the central region and the leads in the derivation 
of the current formula. Given a general basis set, {&}, 
(which we refer to as the direct basis) for the finite dimen- 
sional Hilbert space, H, there exists a dual basis, 
with the property 

(4>i\ct>j) = S ij . (18) 

The vectors of the dual basis can be represented explicitly 
in terms of the direct basis, 

We shall make the general convention that indices 
marked with a bar refer to the dual basis. From the 
expansion l|19fl it follows that the overlap matrix of the 
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dual basis is simply the inverse of S, i.e. 



s J j=(<k\<t>j) = (s- 1 h- 



Before coupling: 



(20) 



Here it is to make a connection with the work of Fransson 
et al., who introduce a set of creation/annihilation oper- 
ators fulfilling the anti-commutation relations {cj,ct} = 
(5 _1 )i 3 (this should be contrasted with Eq. 0). How- 
ever, as can be seen from Eq. (1) of Ref . ITsI these are sim- 
ply the creation/annihilation operators of the dual basis, 
and thus the formalism of Fransson et al. is consistent 
with the one presented here. 

As a final observation we note that the following rela- 
tion holds for any of the single-particle GFs 



(21) 



That is, the GF matrix in the dual basis equals the over- 
lap GF in the direct basis. 



III. CURRENT FORMULA 
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FIG. 1: Schematic diagram of the system used to derive the 
current formula. Before the coupling between the three re- 
gions is established, the leads are in equilibrium with chemical 
potentials \ll and fiR, respectively. Upon coupling a current, 
/, starts to flow as the system approaches a common equilib- 
rium. Interactions are allowed inside the central region (C). 



In this section we derive a formula for the electron 
current through a spatially confined region possibly con- 
taining interactions. The derivation follows the original 
work of Meir and Wingreeni^, but is here extended to 
the case of a nonorthogonal basis set. Atomic units will 
be used throughout. 

We consider the transport of electrons through a sys- 
tem which can be divided into three regions (see Fig.[T]l: 
A left lead (L), a right lead (R), and a central region (C). 
For times t < to the three regions are uncoupled, and the 
leads are in thermal equilibrium with chemical potentials 
Hl and /i^, respectively. When the leads, which we as- 
sume to be macroscopic yet finite in size, are coupled to 
the central region a current will start to flow as the sys- 
tem approaches a common equilibrium. The macroscopic 
size of the leads ensures that a steady state with a con- 
stant dc-current will exist for a considerable time before 
the system reaches equilibrium and the current dies out. 
It is the determination of this steady state current that 
we address in the following. 

As a basis of the single-particle Hilbert space, we take 
a set of localized functions each of which can be assigned 
to exactly one of the three regions, e.g. by the center 
position. We denote the elements of this basis by <fi a i, 
where a £ {L, C, R} specifies the region and i enumerates 
the basis function within the region. Since the 4> a i are all 
localized, we can safely assume that there is no overlap 
between two basis functions belonging to different leads. 
Indeed, the central region can always be enlarged by part 
of the leads until this condition is fulfilled. We thus have 
Shi,Rj = and Sm^Lj = for all i,j, which we write 
compactly as Slr = Srl = 0. 

The main difficulty in the derivation of the current 
formula is caused by the nonorthogonality of the basis 
functions in the central region and those in the leads. 
To overcome this problem, we replace the basis functions 



in the central region by the corresponding functions of 
the dual basis, i.e. <f>ci — ► 4^1- The original basis func- 
tions are maintained in the leads. We shall refer to the 
resulting basis, {(4>Li), (<%l)> (</>«)}) simply as the new 
basis. The introduction of the new basis amounts to a 
re-definition of the subspace associated with the central 
region, such that all three subspaces become orthogonal. 
It should be stressed that this is merely a basis change, 
and thus the full Hilbert space as well as any calculated 
physical quantity remains unaffected. The nonorthogo- 
nality of basis functions within a given region presents no 
serious problem, and can be handled by the techniques of 
the previous section. In the following all matrix quanti- 
ties referring to the new basis will be denoted by a tilde. 

The electronic Hamiltonian, H, consists of a non- 
interacting part, h, and a part containing interactions, 
Vint- The physical nature of the interactions is not im- 
portant, however, we assume that Vi n t only affects elec- 
trons located in the central region. In terms of the new 
basis, the matrix associated with the non-interacting part 
of the Hamiltonian has the generic shape: 

_ / h LL h w \ 
h = h m hec hc R ■ (22) 
V h RC h RR J 

Here, for example, the matrix element hci.Lj — 
is given by (in the coordinate representation), 

hci,Lj = J dr0^(r)[-iA r + «(r)]^(r), (23) 

where v(r) is the sum of all external potentials acting on 
the electron system. We shall not be concerned about the 
specific form of v(r), as this has no importance for the 
general treatment presented here. The non-interacting h 
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is further divided into the two terms, ho and h coup : 

H = ho + h coup + Vint, (24) 

where 

ho = ho+ho+h$ (25) 

ae{L,R},i,j i,j 

describes each of the regions separately, and 

h C oup — ^ ^ fya,ajC^7j£aj 4"H.C, (27) 
ae{L,R};Lj 

provides the coupling. As described in Sec. Ill Al the ma- 
trix f) is given by f) = S hS -1 , where S is the overlap 
matrix in the new basis, 

/ S LL \ 

S=\ (S-^cc . (28) 

V o o Srr J 

The form of S within the central region follows from 
Eq. (|20l) . Note that in the new basis there is no overlap 
between basis functions belonging to different regions of 
the system, i.e. the three subspaces are orthogonal. As 
a consequence, all creation/annihilation operators of a 
given region commute with all creation/annihilation op- 
erators of the other regions. This property of the new 
basis is crucial for what follows. 

For times t < to the three regions are decoupled and 
the state of the system is described by the density oper- 
ator 

PO = -j-p%PoPo®Pext, (29) 

where Pq * = e ~^-^) and z = Tr[p° L p c p% ® Pext ], 
while pcxt describes the state of possible external degrees 
of freedom such as phonons or magnetic impurities. No- 
tice that the order of the density operators in Eq. i|29|) 
plays no role since they all commute due to the properties 
of the new basis. 

After the coupling has been switched on at time to, 
the state po will evolve according to the full Hamiltonian 
of Eq. i|24|) . At the later time t\ the system has reached 
equilibrium and the current has died out. We assume 
that there is a time interval [to + At; t\ — At] during which 
the current through the system stays constant, i.e. the 
system is in a steady state. In the following we consider 
times t for which the system is in the steady state. 

The particle current from the left region into the cen- 
tral region at time t is given by the time derivative of the 
particle number in lead L^i 

I L {t) = ^(N L {t)) = i{[H,N L ][t)). (30) 



The second quantized form of Nl follows from Eqs. (@J 
and JSJ) when is the orthogonal projection onto the 
subspace spanned by the basis functions of lead L, 

#i = X>^)« c W (31) 

Note that only basis functions of the left lead occur in the 
expression for Nl- If we had used the original basis, the 
expression would also contain creation/annihilation oper- 
ators of the regions C and R. Since the interaction, Vi n t, 
by assumption contains creation/annihilation operators 
of the central region only, the commutator in Eq. I|3U|I 
vanishes for all terms in H except those coupling L and 
C, i.e. 

= J Tt[1) L cG< l (lo) - G< c (u)i, C L]duo (32) 

where the commutation relations Eqs. ill I- ill have been 
used and the lesser GF defined in Eq. has been in- 
troduced. In the second line we have assumed that in 
the steady state the GFs depend only on the time dif- 
ference t — t', and moreover that the steady state exists 
for sufficiently long that boundary effects associated with 
the switching on of the coupling and levelling out of the 
current can be neglected when performing the Fourier 
transform. These conditions can always be fulfilled by 
increasing the size of the leads. In the following explicit 
reference to the w-dependence will sometimes be omitted 
to simplify the notation. 

The lesser GF can be obtained from its contour- 
ordered counter parts via analytic continuation as de- 
scribed in Ref. 21. Treating h COU p and V%nt perturba- 
tively, the rules for perturbation theory in a nonorthogo- 
nal basis (see Sec. lIIB(l lead to the following Dyson equa- 
tions for the lesser GF matrix 

G CL = G CcbcLg°LL + G CC^CLg° L t (33) 

g lc = 9ll^lcG^c +9ll^lcGcc, ( 34 ) 

where g\ L is the GF of the uncoupled left lead, i.e. the 
GF defined by ho and po- Here it is important to real- 
ize that this is not an equilibrium GF, since po involves 
different chemical potentials and therefore is not an equi- 
librium state. However, since the CLi,c Li commute with 
the corresponding operators of region C and R (due to 
the properties of the new basis), g\ L ^ s m ^ ac ^ equal to 
the GF defined by and pg . Since the latter describes 
a system in equilibrium the fluctuation-dissipation theo- 
rem provides the relation 

9ll M = ~h M \9ll H - 9ll Ml > (35) 

where /l (uj) is the Fermi distribution function of the left 
lead. The relation l|35|l introduces the equilibrium dis- 
tribution of the lead into the current formula. Since the 
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strict validity of this relation, as explained above, relies 
on the orthogonality of the three regions, we again see 
the importance of working in new basis. 

For the contour-ordered GF matrix of the central re- 
gion we have the Dyson equation 

Gcc — Gq C + Gq C [El + £7? + T,int]Gcc, (36) 

where and £ a are self-energies due to the interac- 
tions and the coupling to lead a, respectively. Note that 
the former also contains contributions from the coupling 
since a complete separation of the diagrams related to the 
two perturbations is not possible. The contour-ordered 
self-energy matrix due to the coupling to lead a is given 
by 

Z a (tu) = t) Ca g° aa (uj)t) a c- (37) 

It is useful to introduce the coupling strength due to lead 
a, 

r a H=*K(w)-E£(u;)]. (38) 

Inserting G CL and G^ c from Eqs. l|3~3*j) . (|3~4*|l into the ex- 
pression for the current, Eq. (|32|> . and symmetrizing, 
/ = (II — I R )/2, we arrive, after some algebra, to the 
desired current formula 

I = \j Tr[(f L -f R )G< c 

+ (f L (u)T L - f R (u)f R ){G cc " G a cc )]du. (39) 

Eq. (|39[l is formally equivalent to the corresponding for- 
mula valid in an orthonormal basis^, and indeed it re- 
duces to it when the basis is orthonormal. However, 
it should be remembered that all quantities entering 
Eq. (|39|) refer to the new basis (as indicated by the tildes), 
and that the coupling strengths, T l and T Rl involve f) in- 
stead of h. 

As the dual basis in most cases is not explicitly known 
it is desirable to re-express Eq. (|39|) in terms of the origi- 
nal basis, that is, in terms of untilted quantities. For the 
GFs of the central region we have the simple relation 

G cc (cu) = Gf^) = GZ c (u), (40) 

where the second equality follows directly from Eq. (|21[) . 
As for the Ts we note that this relation holds in particular 
when Vint = 0. In this case Eq. (|16[) establishes that 
®cc = ~ h)~ x ]cc- On the other hand G r cc can 

be obtained from the Dyson equation (|36fl using that 
G^qq = [oj + S^ — S^hccS^Q^ 1 which in turn follows 
from Eq. 114fl . In terms of self-energies we thus have 

<3 r cc = [uj+Scc-hcc-Zi-^y 1 (41) 

G r C c = [ w+ <Scc ~ ScchccScc ~ ~ s r] 1 ( 42 ) 
where 24 

= {u + S Ca - h Ca )[u + S aa - h aa }~ 1 (uj + S a c - h aC )- 

(43) 



By equating Eqs. (|41(l and (|42|l . it follows after some 
matrix algebra, that 

-> r Q as 77 -> 0, (44) 

where r Q = — S^] and ir\ is the imaginary part of 
uj + . Eq. 144|) holds, of course, also when interactions are 
present. 

With equations l|4(J|l and (|44|) we can state our main 
result, namely an expression for the current in terms of 
quantities which are all evaluated in terms of the original 
nonorthogonal basis: 

I = \J Tt[(T l -T r )(3< c 

+ (f L (oj)T L - f R (cu)T R )(® cc - <S a cc )]dtu. (45) 

A. Non-interacting electrons 

In the special case of non-interacting electrons the 
Keldysh equation^ for the lesser GF of the central re- 
gion in combination with Eqs. (|35|) . l|37l) . and l|38l) yields 

G cc = G uc &< + X<)G a cc (46) 
= G uc [f L {u)T L + f R (cj)T R ]G cc . (47) 

(The Keldysh equation relies only on the temporal prop- 
erties of the various GFs and is therefore valid in any 
single-particle basis). Substituting this relation into 
Eq. (H| and using the general result G r -G a = G>-G<, 
we obtain the following Landauer-type formula for the 
current, 

I = J [f L {uj) - f R (L,)]Tr[G r cc T L G a cc r R ]dcj. (48) 

By virtue of the identities l|4l)|l and (|44|l we can again 
re-express Eq. I)48|) in terms of quantities of the original 
basis: 

I = J \f L (u) - fR(u)]Tr[G r cc r L G% c r R ]du. (49) 

This is the celebrated "trace-formula" which has been 
widely used for numerical calculations of coherent trans- 
port. The same formula has previously been derived by 
Xue et al. using a somewhat different approach involving 
transformations to a real space representation. 

B. Interactions 

In the presence of interactions in the central region, 
the current formulas l|39|) , H45[) are exact provided the full 
interacting Green's function, Gcc {=®cc), is known. In 
this section we address the evaluation of the interacting 
GF within perturbation theory given a nonorthogonal 



7 



basis. To avoid any confusion we denote the GFs eval- 
uated in the presence of the coupling to the leads, but 
without the interaction, by the superscript "ni". In the 
following we use the idea of the "new basis" introduced 
in the beginning of Sec. IIIII without further explanations. 

1. Electron- electron interactions 

Assume that the electrons located in the central re- 
gion can interact through a two-body potential, V^ 2 \ 
The direct basis of TiS 2 ^ consists of the tensor products 
{4> a i ® 4>/3j}j where a, (3 € {L, C, R}, while the new ba- 
sis is obtained by replacing (ftci by its dual <^j. The 
limitation of interactions to the central region means 
mathematically that the matrix element Vdpj^kSi = 
{4>ai ® (frpjlV&'l^k ® <j>si) is non-zero only when all 
a,/3,j,S = S. Since the single-particle basis functions 
used in the lead regions are the same in the new and the 
original basis, the same holds for the matrix V, as well as 
for the matrix 2} = {S^^ViS^)- 1 . The first claim 
can be verified by expanding <j>^ in terms of the direct 
basis (<f>ci), (4>Ri)} according to Eq. ljT§|) . and the 
second then follows from direct calculation using Eqs. (J7J 
and ||2EJ). 

As described in Sec. Ill Al the second quantized version 
of i n the new basis reads: 

Vint = ^2^CiCj,ckcic^-cL.c^ I c m . (50) 

ijkl 

By direct calculation it can be shown that 93 = V, and 
thus 

Vim = ^2 V ^C3,CkClC^-Cy-C mCm . (51) 
ijkl 

Treating Vi n t as a perturbation to the non-interacting 
Hamiltonian, h = ho + h coup , the corrections (diagrams) 
to the non-interacting GF of the central region should 
thus be evaluated according to the usual Feynman rules 
using V as the coupling matrix an d G™ c as the free prop- 
agator, see last paragraph of Sec. Ill Bl From Eq. I|40|l we 

get that G™ c = ^c^C anc ^ consequently we have estab- 
lished the perturbation series for the central object of 
the current formula, Gcc, m terms of V and <3™ c , i.e. 
quantities referring to the original basis. 

2. Electron-phonon interactions 

The interaction between the electron system and a sin- 
gle vibrational mode is, to first order in the ion displace- 
ments, described by an operator of the generic form: 

V mt =v^{v)[b + b\ (52) 



Here W , b are creation and annihilation operators acting 
on the phonon system, and v^ 1 ' (r) is a one-body poten- 
tial acting on the electrons, i>W is obtained by differ- 
entiating the electron-ion potential, V e i-i on (r, R), with 
respect to the ion-coordinates in the direction of the vi- 
bration under consideration, see e.g. Ref. |22|. The re- 
striction of interactions to the central region implies that 
v ai,i3j = {4>ai \4>l3j) ls non-zero only when a, (3 = C. 
Physically this means that the vibration does not distort 
the potential felt by an electron outside the central re- 
gion. According to Eq. 10, the matrix that is relevant 
for the second quantized form of v^- 1 ' in the new basis, is 
D = S~ l vS~ l . Using the expansion (|19f) and the block 
diagonal form of S (see Eq. I|28ll ). it is straightforward 
to establish that = v. The second quantized version of 
Eq. l(52"|) thus reads 

V int = J2 v Ci,Cjcy>cj[b + b^. (53) 

ij 

Considering Vi n t as a perturbation to the non-interacting 
Hamiltonian, h = Iiq + h coup , the diagrams generated by 
the interaction should thus be evaluated using v as the 

coupling matrix and G™c ( = ®cc ) as the f ree propaga- 
tor, see last paragraph of Sec. Ill Bl 



IV. SUMMARY 

The use of localized basis sets in electronic struc- 
ture calculations calls for general formulations of applied 
physical theories, taking the nonorthogonality of the ba- 
sis into account. In this paper, I have presented a general 
form of the second quantization formalism and Green's 
function theory which is valid in a nonorthogonal basis, 
and used it to obtain a nonorthogonal version of Eqs. (6) 
and (7) in Ref. for the current through an interact- 
ing electron region. The main problem in deriving the 
generalized current formula, namely that the lead sub- 
spaces are not orthogonal to the central region subspace, 
was solved by replacing the basis functions of the central 
region by the corresponding elements of the dual basis. 
This simply amounts to a basis change, under which the 
central region subspace becomes orthogonal to the leads. 
The nonorthogonality of the basis functions within each 
of the regions was handled by applying the generalized 
form of the second quantization and Green's function for- 
malisms. Finally, the appropriate nonorthogonal form of 
the perturbation expansion for the Green's function was 
established for the case of electron-electron and electron- 
phonon interactions in the central region. 
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